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Abstract 

Using an effective potential approach, we compute two-loop radiative corrections to 
the MSSM lightest CP-even Higgs boson mass M^o to 0{a^) for arbitrary left-right 
top-squark mixing and tan/5. We find that these corrections can increase M^o by as 
much as 5 GeV; assuming a SUSY scale of 1 TeV, the upper bound on the Higgs 
boson mass is M^o ^ 129 ± 5 GeV for the top quark pole mass 175 ± 5 GeV. We 
also derive an analytical approximation formula for M^o which is good to a precision 
of < 0.5 GeV for most of the parameter space and suitable to be further improved 
by including renormalization group resummation of leading and next-to-leading order 
logarithmic terms. Our final compact formula admits a clear physical interpretation: 
radiative corrections up to the two-loop level can be well approximated by a one-loop 
expression with parameters evaluated at the appropriate scales, plus a smaller finite 
two-loop threshold correction term. 
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1 Introduction 



It is difficult to overestimate the importance of the experimental discovery of the Higgs boson. It 
would not only help us to elucidate the dynamics responsible for electroweak symmetry breaking 
but it will most probably offer also an important clue as to the nature of the Physics beyond the 
Standard Model (SM). 

The paradigm for this new Physics, granted that a fundamental scalar drives electroweak sym- 
metry breaking, is the Minimal Supersymmetric Standard Model (MSSM) the most econom- 
ical extension of the Standard Model that incorporates (softly-broken) Supersymmetry (SUSY). 
In spite of the uncertainties related to the origin of supersymmetry breaking (and therefore of the 
masses of the so far undetected supersymmetric particles), it is well known that the MSSM pre- 
dicts the existence of a hght Higgs particle with mass below about 135 GeV (this bound depends 
sensitively on the top quark mass one uses; our present calculation intends to set a precise and 
firm bound). Unlike the case of the Standard Model (in which the mass of the Higgs boson is an 
unknown parameter), the mass of the light CP-even Higgs boson of the MSSM is calculable as a 
function of other masses of the model. A precise calculation of that mass is of prime importance 
for Higgs searches at LEP, Tevatron and the LHC, and is the topic of this paper. 

We recall at this point that the Higgs sector of the MSSM consists of two SU (2) doublets. Hi 
(which gives mass to down-type quarks and charged leptons) and H2 (which gives mass to up-type 
quarks). The vacuum expectation values (^1,2) of these doublets break the electroweak symmetry, 
after which, the Higgs spectrum contains two CP-even scalars {h^ and H^; with rriho < itlho), one 
CP-odd pseudoscalar (A°) and a pair of charged Higgses {H"^). At tree-level, the masses, couplings 
and mixing angles of these particles are determined by one unknown mass parameter (say m^o) 
and the parameter f3, which measures the ratio V2/vi{= tan/?). In the limit m^o ^ Mz all the 
Higgs particles except have masses ~ m^o and rearrange in a complete SU{2) doublet almost 
decoupled from electroweak symmetry breaking, while /i" remains light with m^o < Af| cos^ 2/3 
and has SM properties. This bound (which applies for any value of m^o and is saturated for 
tuao ^ Mz) is extremely important: it represents the limit that experimental bounds should 
reach to falsify the MSSM. In fact, the present experimental bound from LEP, including the 
latest data with up to ^/s = 202 GeV, is rriho > 107.7 GeV (for large m^o, case in which the SM 
limit is applicable; the limit falls to ~ 91 GeV for smaller m^o), which is well above this bound. 
This is not yet conclusive evidence against the MSSM because it does not take into account the 
radiatively corrected form of the mass bound. 

Radiative corrections to m\o have been computed using three different techniques (or combina- 
tions of them): effective potential method P, |, |^, ||, 0, direct diagrammatic calculation P, |], 
and effective theory (or renormalization group) approach ||^, |^ |TT|, |T^, |T^. The full one-loop 



radiative corrections to m^o have been computed diagrammatically. The most important of these 
corrections come from top quark/squark loops and are given by 

27r^t;^ mf 

where rrit is the top quark mass, an average top-squark mass and v"^ = vf + V2 = (246 GeV)^. 
This correction can be very large if mj ^ rrit, and in such case rrihO can evade easily the current 
experimental lower bound. This important 0{at) logarithmic correction to the dimensionless 
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ratio m1o/mf [here at = hf/{47c), where ht is the top-quark Yukawa couphng] can be most easily 
reproduced using renormahzation group (RG) techniques. In addition, there is a finite (non- 
logarithmic) correction which may also be important, and which depends on the details of the 
top-squark spectrum. This correction is (assuming again for simplicity degenerate soft masses for 
the top-squarks) 

where Xt = At + n cot /? is the top-squark mixing parameter. At the soft trilinear coupling as- 
sociated to the top- Yukawa term in the superpotential and /i the supersymmetric Higgs mass 
parameter. Correction (^ is maximized for X'f = 6m? (the so-called 'maximal- mixing' case). 
When using one-loop equations like (|lD and (H) to compute the Higgs mass one has to decide 
whether to use on-shell (OS) or running values for the mass parameters that enter such formulae 
(and if the latter, at which scale to evaluate them). The difference between two such choices 
is of higher order, but can be non-negligible, especially because of the m^-dependence of Am^o- 
Although RG techniques can be used to make an educated guess of the scale at which those mass 
parameters should be evaluated (see e.g. [T^), a precise answer to such questions could only be 
unambiguously given by a two-loop calculation like the one we perform in this paper. 

At two loops, radiative corrections to m^o depend not only on the large top- Yukawa coupling 
but also on the QCD coupling gs. It is reasonable to expect that the dominant two-loop corrections 
will be of order (9(asat[ln(m?/m^)]'^) and 0{a'f[ln{'m'^/m'f)]^), k = 0,1,2. Terms with k = 2 are 
the two-loop leading logarithmic contributions and can be obtained by RG techniques using one- 
loop RG equations; no true two-loop calculation is required and RG resummation will take into 
account such leading-logarithmic (LL) corrections to all loops. The k = 1 terms are the two- 
loop next-to-leading-logarithmic (NTLL) corrections, which can be obtained (and resummed to 
all loops) with two-loop RG equations. Finally, the two-loop non-logarithmic terms {k = 0) can 
be interpreted in the effective theory language as threshold corrections (at the supersymmetric 
scale set by the mass of the top-squarks) and require a genuine two-loop calculation; they simply 
cannot be obtained from RG arguments. 

The status of these higher-loop calculations of the radiatively corrected m^o is the following. 
Higher-order logarithmic corrections were included in studies which used RG techniques almost 
since the dramatic impact of radiative corrections on ruho was first recognized. Hempfiing and 
Hoang were the first to perform a genuine two-loop calculation of m/jo which also included 
non-logarithmic terms. They computed the dominant two-loop radiative corrections [to 0{as<yt) 
and 0{a^)] in the case tan/3 ^ 1 and zero top-squark mixing. Their computation also included 
the most important logarithmic corrections, which could be alternatively incorporated by RG 
resummation from one-loop results, as done e.g. in Ref. In this last paper it was also pointed 
out that by a judicious choice of the renormahzation scale at which to evaluate one-loop corrections, 
the higher order logarithmic corrections could be automatically taken into account. A similar idea 
was later implemented in [^, |T^ to write down simple analytical approximations for the radiatively 
corrected m^o, obtained by iterative integration of RG equations. 

Besides being limited to a particularly simple value of tan (3, the calculation in Ref. [^ missed 
the sizable impact of non-zero top-squark mixing in two-loop effects, that is, higher order correc- 
tions to the contribution written down in Eq. (Q). Such corrections were first included to order 
0{as(yt) in the diagrammatic calculation [Q, and by the effective potential method in Ref. [0. 
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The effect of these corrections is to shift the values of Xt that give maximal mixing, change the 
corresponding Higgs mass by up to ~ —10 GeVQ and introduce an asymmetry in the dependence 
of ruho with the sign of Xt. This two- loop top-squark mixing dependent correction was also ex- 
phcitly isolated recently by the present authors in Ref . , which uses effective potential plus RG 
techniques. Besides confirming the independent diagrammatic results of Ref. |^ we clarified the 
relation of these calculations to previous ones (in particular matching results expressed in different 
renormalization schemes; see also 13). We also derived a compact formula for the Higgs mass 



(in the spirit of |T2|, 0) which took into account the most important radiative corrections, and 
used RG techniques to include in a compact way two-loop LL and NTLL corrections. With the 
0{as<yt) radiative corrections organized in this way, we find that the mixing-dependent genuine 
two- loop threshold corrections are generally small 3 GeV). 

Nevertheless, the large computing effort just reviewed did not exhaust the potentially impor- 
tant radiative corrections: the two-loop 0(aj) top-squark-mixing-dependent corrections to rUhO 
remained unknown to this day, while it is clear that they could compete in principle with the 
0{as(yt) contributions. The purpose of this paper is to complete the calculation performed in 
[^, ^ by using effective potential techniques (plus RG techniques) to compute such C(a^) con- 
tributions for general top-squark mixing parameters and any value of tan/3. The results in this 
paper can be considered the most complete and accurate approximation to rriho presented in the 
literature. 

The structure of the paper is the following: the next Section describes the strategy of our 
calculation and presents some analytical formulae for m/jO, obtained in the limit of ^ nit. 
Section 3 goes one step ahead implementing the RG-improvement of such approximations and, 
in doing so, clarifies the organization of the higher order radiative corrections. This procedure 
is not only important to provide a clearer physical picture in connection with the effective field 
theory but also to classify those corrections calculated in Sec. 2 into a numerically dominant and 
compact part plus smaller finite threshold correction terms. In Section 4 we present our numerical 
results for the Higgs mass, illustrate the size of the new corrections and check the validity of our 
analytical approximation formulae. We draw some conclusions in Section 5. 

Several appendices are devoted to technical details of different aspects of the calculation. 
Appendix D is worth special mention as it contains the two-loop 0{a^) MSSM effective potential 
used as starting point of our calculation and first computed in this paper. 



2 CP-even Higgs boson masses to two-loop order 



The momentum-dependent mass-squared matrix for the CP-even Higgs bosons of the MSSM in 
the interaction eigenstate basis hi, /i2 is 



mis} + m\,cl + AMisi/) 



(3) 



where S/? = sin/3 and cp = cos/3. The mass parameters niz and uiao are the (scale-dependent) 
running masses of the Z-boson and CP-odd Higgs boson they are related to the on-shell masses 



^This correction is relative to the one-loop mass using on-shell parameters. The size of the correction would be 
much smaller if running parameters are used in the one-loop formula (|l|). 
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Mz and M^o (we use capital letters to distinguish on-sliell parameters from running ones) by 



rrir 



M| + RenL(M; 



^1) 



m 



AO 



Re n^A(Mlo) - s 



Vl 



- c, 



• T2 

'V2 



(4) 



where is the transverse part of the Z-boson self-energy and 11^^ the A'^-boson self-energy, 
Ti, T2 are the tadpoles of the CP-even (real) fields hi, h2- Their explicit one-loop expressions can 
be found e.g. in Ref. ||10 



In (|^), AAi^ stands for the contributions from radiative corrections. They are 



^mUp 



(5) 



where Uij is the self-energy matrix of the Higgs fields hi and h2. The masses, rriho, m^^o, of the 
two CP-even Higgs bosons are then obtained from the real part of the poles of the propagator 
matrix, 



Det 



M 



hy^ho^Ho] 



. 



(6) 



The radiatively corrected mixing angle a is obtained as the angle of that rotation which diagonal- 
izes M.\ (for some choice of p^, say 



tan2Q; 



(MDii - (Ml 



(7) 
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Computing Higgs boson masses to a certain order of perturbation theory then requires calculating 
the self-energies and tadpoles in Eqs. @ and to that order. 

In the effective potential approach |^, ^, |, ^ |^, self-energies and tadpoles can be calculated 
as derivatives of the Higgs potential V according to: 



T, 



dV{hi,h2) 
dhi 



n.,(o) 



hi=vi,h2=V2 



d^V{hi,h2 



dhi dhi 



(8) 



hi=vi,h2=V2 



Note that the H's obtained from derivatives of V have zero external momentum. 

In the limit M^o 3> Mz, the lightest CP-even Higgs state lies along the direction of the breaking 
in field space [0, that is, a —>■ P — 7i/2 + 0{m%/m\), and its radiatively corrected mass has a 
very simple expression 



M? 



4mf 



d 



hO 



dmf 



V-Re Uhhimlo) + Re H;,;,(0) , 



(9) 



which is exact up to corrections of order (9(m|/m^o)-Q 

In Eq. (P), V is the projection of V{hi,h2) along the light Higgs h = hiCp + h2Si3: V{h) = 
V{hi hc/3, h2 — > hs/s). Then V{h) can be expressed as a function of rrit using h — > mt\/2/ {htSp). 

^ This formula can be proved as follows: If m^o S> to^j ^ ~* /3— 7''/2, and we can therefore use the approximatfon 
Am^o ~ AMliC^ + AA^Ij*^ + SAAl^j'^/SC/j, up to higher order terms in m^/m^o. Observing that the potential 
V depends on the fields hi and /12 only through (field-dependent) top quark mass and the off-diagonal elements of 
the top-squark mass-squared matrix and using (|^), we can easily express the partial derivatives of V in terms of 
the total derivative in (^. A similar formula was already used in 
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We decompose V in its n^'^-loop pieces Vn (explicitly given in Appendix D) as V = Vq + Vi + V2. 
The tree-level part Vq is the only one in which we keep non-zero electroweak gauge couplings. We 
approximate the one-loop part Vi by its 0{at) piece coming from top quark/squark loops. The 
two-loop part V2 is approximated by V2S + V2t, where V2S is the 0{asat) part and V2t the 0{a^) 
one. 

Next, Uhhip^) is the light Higgs self-energy at external momentum p, related to the self-energies 
of hi^2 by 

Uhhip^) = Uu{p')sl + Il22{p')cl - 2Uu{p')so.c^ . (10) 

Notice that the self-energy difference in (^ involves non-zero external momentum and would 
require a diagrammatic two-loop calculation. However, throughout this paper we work in the 
approximation of neglecting in the radiative corrections all couplings other than hf or g^. In 
that case, realizing that at tree level rriho depends only on electroweak gauge couplings while its 
dependence on ht appears only at one-loop, we can write 



m 



_d_ 

dp"" 



which gives rise to 0{af) contributions from the one-loop 0(at) self-energy Hhh- 

In Section 4 we present the numerical results of such procedure for the two-loop potential with 
0{asat) and C(a^) corrections included. The general expression for the 0{asat) potential was 
first computed in while the complete 0{a^) terms were still missing. Both contributions to V 
are given in Appendix D. 

It is useful, both for a better understanding of the numerical results and for practical applica- 
tions, to derive an analytical expression for the light Higgs mass in the case of a large hierarchy 
between the supersymmetric scale and the electroweak scale (say when the SUSY scale is of order 
1 TeV). Such limit is interesting because it maximizes the radiative corrections to m^o (so that 
it corresponds to the most pessimistic scenario for Higgs searches; the case one should be able 
to discard to rule out the MSSM), and at the same time simplifies the structure of the radiative 
corrections, avoiding the proliferation of a multitude of different supersymmetric thresholds. 

We consequently assume now that all supersymmetric particles have roughly the same mass 
Ms ^ Mz- In more detail, focusing on the particles relevant for the radiative corrections to 



m/jo, we take equal soft masses M 



Q 



u 



Ms for the top-squarks (with diagonal masses 



M| + m'f). The two eigenvalues and mixing angle of the top-squark squared-mass matrix 



are then 



■rrir 



mr + mtXt 



t , 



m 



t2 



rrir 



mtXt 



(12) 



We also take the same mass Ms for the gluino and the pseudoscalar Higgs [this means in particular 
that we can use Eq. (|]) for the light Higgs boson]. In principle we admit the possibility that 
the /i parameter could be smaller than Ms-, in which case we expect that one chargino and two 
neutralinos will have masses ~ below the common supersymmetric threshold. In this situation, 
which broadly corresponds to the case of a common heavy SUSY scale, we find that, using the 
operator 
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the different parts entering are 
VlV^ = m| 008^2/?, 



vlv. 



2t 



3m^ 



In ^ + * 



m 



In 



2^1 



2 In 



Ml 



12 

2^1 



+ 21n 



2^ 



+ In 



1+ -l + 21n— | + 21n 



3 J 12 



(14) 
(15) 

(16) 



9 In^ ^ - 6 In In ^ - 3 In^ ^ + 2[3/2(/i) - 3/i (/i) - 8] In ^ 



Ml 



efiHi- In 



M| 



2(4 + /i2)/i(/i)+4/3(/i)- 



TT 



+ 



+ 



(33 + 6/i^) In - 10 - 6/}^ - 4/2(/i) + (4 - 6/i=^)/i(^) 
-4(7 + /i^) In + 23 + 4/i2 + 2/2(/i) - 2(1 - 2ix^)Mfi) 



Ml 



Ml 



31n^^ + 71n^ - 41n^ -3 + 60K + 



47r2 



Mi 



+ 12-24ir- 181n— f X/- [3 + IQK ~ 3\n-^ ]{AXtYt + 



Ml 



+ 



11. Ml 



Ml 



6 + — In-f X,^+ 4 + 16ir-21n— f Xj^Fj 



14 



Ml 



-= + 24i^ - 3 In ^ ) XfY^ - I ^ + SiT - - In ^ ) X'^Y^ 



19 



1, mV 



-4V>2 



12 



(17) 



The notations used are Xt = At + (1 cot /5, = — /i tan /5, with reduced parameters z = z/Ms, 
and (see Appendix A) ~ —0.1953256. We also use the following non-singular functions of fi 



/2(/i) 
/3(/i) 



l-/i2 
1 

l-/i2 



In/i^ 



1 + 



1 — jj,' 



■ In /t^ 



-1 + 2/i2 + 2/i^ 



2^2 



In/i^ ln(l - fi^) + Li2{fi^ 



7r 
"6 



— fi^ In /t^ 



(1^ 



with /i(0) = 0, /2(0) = 1, /3(0) = -7rV6 and /i(l) = -1, /2(1) = 1/2, /3(1) = -9/4. 

Finally, the correction for non-zero external momentum in Eq. (P) is given by (see Appendix C) 



Re 



^m,o.,(31n- + 2-- 



X? 



(19) 
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The parameters that appear in these expressions are running parameters, evaluated in the 
DR -scheme and satisfy MSSM RG equations. In fact, it can be checked that the (physical) Higgs 
mass, given by Eq. (|^), is renormalization-scale independent (up to two- loop order), as it should. 
This scale independence is at the root of the RG-resummation procedure discussed in the next 
section. It is evident that for different values of the renormalization scale, the magnitude of the 
two-loop corrections will change, so that it should be possible to choose the scale in such a way 
that the bulk of the corrections is transferred to the one-loop terms (which depend on the scale 
implicitly). 

Therefore, the magnitude and relevance of the two-loop corrections depends on the definition 
of the mass parameters that enter the one-loop corrections. It is in this respect convenient to write 
down the two-loop expressions just obtained in the particular case in which all mass parameters in 
the one-loop correction are the OS ones. This is also useful to compare with explicit diagrammatic 
calculations. The relationships between running and OS parameters are listed in Appendix C. 
Using them, we obtain for the two-loop correction to mj^o'- 



hO 



Ml 
ml 
Ml 



, f <! -3 In^ ^ - 6 In ^ + 6Xt - 3 In ^ - ^Xf 



Ml 



Ml 



mi 



mt 



7 



3 
4' 



31n2 -f + 131n-f )4 - - - y - 3/i2 - (11 - /i^ + 3fi')j\{fi) 



13 47r2 



+ 



-3(1 - fi'Y ln(l - jj') + 3/2(/i) + 4/3(/i) + 4 [60K + - + 



x: 



1 

+ 6 



61 



-26-9/i^ + 3/i(/i) + 3/2(/i) + y4 



Xf + ^Xf 
t 2 * 



+3(1 - /i') (2 - 3/x')/i(/i) - (1 + 3fi^) ln(l - fi^ 



* 6 



(3 _ iQK - nV3){4XtYt + Y^') + [iQK + ^)xfYt 



+ f + 24A^ + 7rV3]x^Y,' - { ^ + gK + ^ )X^Y,' 



7_ 
12 



TT 



2^3 



+ 2Xt- 



Xf 



94M^-(31n^ + ln4)4x3 



(20) 



We emphasize that this expression gives the two-loop corrections when the one-loop contribution 
([T5|) is expressed in terms of OS parameters, that is. 



OS 



3g'M^ 



i-ioop Svr^M^ 



M? 



In^ + ' 



Mi 



xp^V 1 fxP^V 



Mi 



12 \ Mi 



(21) 



Several features of Eq. ( pOD are worth commenting. First, if we restrict Eq. (|20D to tan/3 ^ 1 
and zero At, to compare with the result of Ref. we find the same logarithmic terms. However, 
the 0{at) finite term is different. In particular, that term is sensitive to the value of the parameter 
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II, contrary to what is stated in Ref. 0. Nevertheless, the result quoted for that finite term in 
Ref. is inside the range we would find by varying /t^ from to 1, and the impact of this 
//-dependence on the final Higgs mass is quite small. 

Second, we see that radiative corrections no longer depend on At and fi in the combination Xt 
that appears through the off-diagonal entry of the top-squark mass matrix: besides the explicit 
dependence on the parameter /i already noticed, the quantity Yf also introduces a different com- 
bination of At and /i. This dependence on Yt originates from the H — t — h and H — t — t diagrams 
of Fig. I 

Third, although roughly speaking the top- Yukawa correction has a small pre-factor 3/16 in 
comparison with the QCD correction, this does not guarantee that the new contributions will be 
negligible compared to the QCD one. In fact, we will see that for two-loop top-squark-mixing- 
dependent corrections of (pOD, the top Yukawa contributions have opposite signs as that of the 
QCD corrections and could be as much as 60% of the latter (see Fig. P). In the next Section, we 
will follow RG methods and reorganize these corrections in the effective theory language, with the 
most important corrections of Eq. (pOD reshuffled in a RG-motivated one-loop formula. 



3 Renormalization group resummation 



Before illustrating in Section 4 the impact of the newly computed corrections on the Higgs mass, 
we show in the following how the use of renormalization group techniques |]Tl], |T3| allows us to 
write the previous complicated corrections [see Eq. (^Of)] in a simpler and more transparent way. 



while at the same time it clarifies the connection to the RG programme, which can be used to 
improve the precision of the mass formula by resummation of higher order corrections. 

We already applied this idea in Ref. to the 0{asOit) two-loop corrections. By a convenient 
(and physically well motivated) choice of the scale at which to evaluate running parameters in the 
one-loop mass correction one can absorb large logarithms in Eq. (pOl) . The RG evolution of the 
parameters is given by the corresponding one-loop RG functions listed in Appendix B. 

We use the following equations to relate supersymmetric running parameters at different scales 
[c/. Eqs. (^) and (^)]: 



m?(Q) =m?(Q') 1 + 



167r2 



16 



-9z 



-hi 

2 * 



Xlsl + Y^cl + cl + 2~2i^ 



In 



Q_ 



12 ' 



(22) 



Xt{Q) = Xt{Q') 



167r2 



where we have used Af 



Xt-sj^ + YtC^p, 



16 

-glMs + ?>hl{Xt + Xtsl + Ytcl) 



In 



Q_ 



12 



A't+fi' 



+ Y^c^p and m^^ + /^^ 



(23) 



Notice 



that, to the order we work, it is sufficient to use these one- loop LL approximations to the full RG 
evolution because we are concerned with parameters that appear in a one-loop order term. 

The Standard Model MS top quark mass Wit and the Higgs VEV v are related to the on-shell 
mass Mt and MSSM VEV v by [c/. Eqs. ( |C.2| ) and (|C.1(]| ), from which relevant terms can be easily 
identified] 



(Q) = M,^ 



1 




(24) 
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v\Q) = v\Q) 



2„2 



(25) 



We also use one-loop LL solutions of the SM RG equations to relate these parameters at different 
scales: 



1 + 



1 



167r2 



12 



v'\Q)=v\Q" 



1 + 



-his^, ] In 
In 



Q_ 



(26) 



(27) 



167r2 g2 

Using the above equations, we find the following compact expression for the Higgs boson mass 
which is one of the main results of this paper 



M;^o = M| cos^ 2(3 + 



The one-loop threshold correction is 



^th "^hO 



3m^(Qth) 

and the two-loop threshold correction reads 



12 



th '"■/lO 



^th "^/lO 



+ 



+ C 



^3^2 

IGvT'^f 2 
2 



* 3 * 12 * 6 *. 
i?o(A) + i?2(/i)X,2 + i?4(A)^' - 1:4X1 



2 i>6 



q Att^ , ^ ^ 

60ir - - + — - (3 + 16ir)(4Xiy; + r/) + (15 - 24ir)X, 



25 



^Xf + (4 + 16K)Xi'i; + + 24K]x^Y,' - ( S + )^'^* 



V 3 



19 



12 



We have used the short-hand notation 

9 vr^ 



+ Gfi' - (11 + 2/i^)/i(/i) + 3/2(/i) + 4/3(/i), 



i?o(/i) = 

/?2(/i) = -ll-/i'[6 + 6/i(/i) + 10/2(/i)] 
R,{p.) = 6 + + Mil) + f2ifi)] . 



The scales required in (^ 



are 



-1/3, 



(28) 



(29) 



(30) 



(31) 



(32) 



Qi = Qth = mi , 

It is a non-trivial check of our calculation that the values of the scales (|3^) required to re-absorb 
the large ln(M|/m^) logarithms in the two-loop corrections are consistent with the ones obtained 
in for the QCD corrections alone. We see, in particular, that the uncertainty found there in 
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the definition of tfie scales Q[ and Qi is liere resolved by the need of absorbing the new radiative 
corrections. 

We still find a somewhat complicated expression for the threshold correction A|^'*m^o, due to 
the fact that we have kept free the /i parameter. Expressions for the two limiting cases of heavy 
/X (/X ~ Ms) and light (/i <C Ms) can be readily derived. In both cases, the resulting threshold 
correction is much simpler than the general case (^0[) and contains no more logarithms. Explicitly, 
for /i <^ Ms we find 

i?o(0) = i?2(0) = -11, i?4(0) = 6, (33) 



and for /i ~ M5: 



2' 



7^2 13 

i?o(l) = 7-y, i?2(l) = -16, i?4(l) = y. (34) 



It is perhaps convenient to make more explicit the connection between our results and those 
obtained in the RG approach (see, e.g., Ref. Il3|). To be concrete, let us assume = Ms so 



that all supersymmetric particles (including charginos and neutralinos) have masses of order M5; 
below that scale, the effective theory is the SM. The light Higgs quartic coupling A at Ms consists 
of a tree-level part plus higher-order threshold corrections which arise from the heavy decoupling 
supersymmetric particles, it can be evolved down to the electroweak scale, say Q = rrit, using the 
SM RGEs; at that scale A is related to the physical Higgs mass. This procedure should reproduce 
all the logarithmic corrections we have found. 

More explicitly, defining j3x = dX/dlnQ'^ , we can write 

X{Qt) = KQt)-r PxdlnQ'. (35) 
J Q=Qt 

We use a special notation for the high and low scales between which we run A to distinguish them 
from other definitions of and that appear in the paper. These quantities are defined by: 

Qt=mtiQt), Qi = mi{Qi), (36) 

i.e., they are the running masses evaluated at a scale equal to the corresponding mass. This is 
the natural definition in the RG approach. 

Making a loop expansion of j3\ in (^) and a further expansion around a particular value of Q 
(say the low energy limit of the running interval, Qi), we obtain to the two-loop order 

A(Qt) ^ A(Q,) - [/3?)(Q0 + /3f (Qi)] In || - ^^^(2*) || + - (37) 

where the one- and two-loop contributions to (3\ are approximated by [neglected all couplings 
other than the strong gauge coupling (73 and the SM top Yukawa coupling gt (= htSjs)] 

tf^ - ^(-^^^ + ^)' 
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(1) 



We note that for a correct two-loop computation it is necessary to retain also the Xg"^ term in j3 
because A gets one-loop contributions proportional to g^. d(3^l^ / dhiQ"^ can be calculated from the 
one-loop RG evolution of gt 



dlnQ2 327r2^""' ^^^^ 
Once \{Qt) is obtained from (|37|) , we extract the physical Higgs mass using the SM relation 
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KQt)v\Qt 



M? 



1 



9t 



(40) 



This correction arises from wave-function renormalization and takes into account the fact that 
the physical mass is defined on-shell, and not at zero external momentum. Its physical content is 
therefore similar to the correction ( pT]) in our effective potential approach. 

According to (|57[), the large LL and NTLL corrections to M^o arise solely from X{Qt). Addi- 
tional radiative contributions in (|40|) , coming from tJ^(Qj) and the wave-function correction factor, 
affect the large logarithmic terms only through multiplication of X{Qt)- It is therefore clear that 
it is enough for our purposes to know these correction factors at one-loop order. Based on this 
observation, we can combine both factors together using (p5|) to write the simpler formula 



(41) 

with Ql = e~^^^Qt, in accordance with (|32D . 

It is now straightforward to show perfect agreement of M^o as obtained from the above expres- 
sion with our results ( ^^ - |30|) . All logarithmic corrections up to two- loops are exactly reproduced 



while the finite part agrees if one uses as boundary condition at the SUSY scale 



with 



4?A 



4 



{gl + gi)cos'2(3 + 6i^X + 6^^>X 



v(2) 



(42) 



3gfiQ 



th 



th J 
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A 



(43) 



To summarize, we find full agreement between our approximate formula (|2^) for the Higgs 
boson mass and the RG-improved mass calculated in the RG (or effective theory) approach, to 
two-loop order. The connection to the effective theory language clarifies the origin of the different 
terms in (pSD , and rewrites them in a very convenient way, absorbing the large (logarithmic) two- 
loop effects in the one-loop correction and leaving behind two-loop threshold corrections which 
are numerically small, as we will see in the next Section. Note that this applies in particular to 
the sizable t op- squark- mixing- dependent corrections of Eq. (pO|) , the bulk of which is transferred 
to the RG- reshuffled one-loop threshold correction of Eq. (|29| ) . 

Knowing the boundary condition, A(Qf), one can integrate (|35| ) numerically by solving a cou- 
pled set of differential equations (describing the two-loop evolution of A, (73, gt), find X{Qt) and 
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use (^T]) to get the Higgs mass. The final result will be the full RG-improved value of Mf-fi and 
will resum LL and NTLL corrections to all loops [numerical integration includes all the terms 
from the expansion around Qt which were neglected in (^7])]. In this respect, note that our com- 
pact formula, Eq. (^), which has been found by requiring that logarithmic contributions are 
correctly reproduced up to two-loops only, contains in fact logarithmic corrections of higher order. 
It can be shown that these higher order logarithmic corrections do not match exactly the correct 
ones (obtained by the RG method) if we use simple one-loop approximations [like those given in 
Eqs. (|2^j25|) ] to evaluate the parameters in (^) at their corresponding scales. However, evaluation 
of those parameters by means of a full numerical integration [similar to that in Eq. (^) for A] 
would correctly take into account the LL (but not the NTLL) terms to all loops. Nevertheless, 
as we will see in the next Section, the error made in neglecting logarithmic corrections of higher 
order is very small for SUSY scales of interest [below Ms ~ C^(l) TeV]. If Ms turns out to be 
significantly larger than that (starting to be in confiict with naturalness criteria), then one should 
revert to the numerical RG integration of A to get a reliable estimate of the Higgs mass. Our 
results for the boundary condition X{Qf) will still be useful in such a case. 



4 Numerical results 

In this section we present numerical results from our two-loop study. For the one-loop analysis we 
closely follow Ref . [^ , which has included complete radiative corrections from the dominant top 
quark/squark sector and the sub-dominant gauge/Higgs boson and neutralino/chargino sectors. 
In what follows, we shall concentrate on two-loop radiative corrections. 

We start by sketching the procedure for this analysis, which is the following: we first take as 
inputs the on-shell mass parameters^ M40, Mt, M~^, M~^ and A^^. From them we can determine 
the values of the corresponding running parameters at any renormalization scale Q. To do this, 
we have to calculate the one-loop self-energy diagrams for Higgses and top-squarks (the latter are 
collected in Appendix C). We also input tan/5 and /i parameters, and convert as{Mz) = 0.118 to 
the MSSM DR running value. Next we calculate in the MSSM the two-loop corrections to the 
CP-even Higgs mass matrix, AJ^fj, from the two- loop potential ( p.5| ) and ( p.6| ) using Eq. (H). 
Numerically, the partial derivatives in these equations are replaced by finite differences in hi, h2, 
i.e. we vary the values of these fields by a finite amount and recalculate the field-dependent 
top-quark mass mf = \h1h\ and top-squark masses m^-^, m^^; mixing angle % from Eq. ( |B.5| ). 
With these new parameters, the two-loop potential is reevaluated and their variations from the 
reference values [calculated at h\ + h2 = (246 GeV)^] are found. Finally, equipped with the 
corrections AAi'fj, we compute the lightest CP-even Higgs boson mass by solving Eq. (|). 

Several approximations have been made to quantities in Eqs. (0-^, in particular we neglect 
all dimensionless couplings other than the top- Yukawa coupling ht and the QCD gauge coupling 
g^. In this way we pick up the dominant radiative effects only, what we term throughout leading 
corrections. We notice that the two-loop self-energy of the Z-boson and the non-zero external 
momentum corrections to two-loop Higgs boson self-energies can be neglected in our calculation 
since all these corrections are higher order effects in the leading approximation. However, we need 
to calculate H^^ to the two-loop level since it has O^agat) and 0{af) corrections and in principle 

■^For the top-squark sector, we can alternatively take as inputs the on-shell top-squark masses and mixing angle. 
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Figure 1: Higgs boson mass M^o vs. the on-shell top-squark mixing parameter X^^. Dotted, 
dot-dashed hnes show one-loop and two-loop 0{asat) results from the program FeynHiggs, cor- 
responding results from our numerical analyses are shown in dashed and solid lines respectively. 




Figure 2: Higgs boson mass M^o vs. the on-shell SUSY scale Ms, for two top-squark mixing 
parameters X^^ = and 2. One-loop mass, two-loop masses to 0{asat) and 0{asat + aj) are 
shown in dashed, dot-dashed and solid lines respectively. 
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Figure 3: Higgs boson mass M^o vs. tan fi for the top-squark mixing parameters X^'^ = and 2. 
Dot-dashed and sohd hnes correspond to two-loop Higgs boson masses to 0{as(yt) and 0{asat+a^) 
respectively for Mt = 175 GeV. Two-loop masses to 0{asat + af) for Mt = 170 and 180 GeV are 
also shown in dotted and dashed lines. 



could contribute to (^) at the same order as AAifj. It is not possible to obtain these self-energies 
in our current approach, and explicit two-loop calculation of the corresponding two-point functions 
are needed. Fortunately, the correction to m/jO from Haa is always numerically negligible for large 
rriAO as can be easily seen from the structure of the Higgs mass matrix (^. That is, (|^) is correct 
for large m^o and we can safely neglect the H^^ corrections. (For m^o ~ niz, a complete two- loop 
calculation of m^o would need H^a-) 

Fig. ID is used as calibration: we compare in it our numerical results for Mho including only 
up to two-loop 0{asat) corrections with the mass obtained by the program FeynHiggs which 
uses the explicit two-loop diagrammatic results of Ref. 0. We choose two sets of parameters^: 
(a) Mao = Mg = M2S = m8^ = Mg = 500 GeV, /i = 200 GeV and (b) M^o = Mg = M2S = 

= Ms = 1 TeV, /i = —500 GeV. For each case, results for two values of tan/5 (1.6 and 20) 
are plotted. We find good agreement (given the fact that they are two independent programs) 
between both one-loop (shown in dotted and dashed lines) and two-loop QCD corrected (shown in 
dot-dashed and solid lines) masses; this shows numerically that the two approaches are equivalent 
to that order. This equivalence is easily understood since the effective potential, as a generating 
functional [|1^ , encompasses all tadpole and self-energy diagrams (as well as all other multi-point 
functions) which are calculated in |^. The effective potential approach is more efficient for the 
purpose of calculating M/jO and much simpler to implement in a Fortran program since it requires 
evaluating only one set of two-loop functions. 

In Fig. we show the Higgs boson mass M^o vs. the (on-shell) SUSY scale Ms, for two values 

^We assume M3 is positive hereafter. For a negative M3 our formulae still apply simply by simultaneous sign 
changes in Xt and l*. 
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Figure 4: Higgs boson masses M/jO vs. X^^. One-loop masses, two-loop masses to 0{asat + af) 
and their approximations are shown in dashed, solid and dot-dashed lines. 
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Figure 5: Higgs boson masses M^o vs. X^^. Two-loop masses to O{as0it) and O{as0it -\- af) 
arc shown in dotted and dashed lines, their corresponding RG-corrected masses are shown in 
dot-dashed and solid lines. 
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Figure 6: Higgs boson masses M^o vs. X^^. Two-loop masses without the top-squark-mixing- 
dependent correction terms of Eq. (pO|) are shown in dotted hues. The corresponding masses 
without the 0{aj) corrections only and the full numerical results are shown in dashed and solid 
lines respectively. 




Figure 7: Higgs boson masses M^o vs. . Two-loop masses without the threshold correction 
term /S.^^m\o are shown in dotted lines. The corresponding masses without the 0{a1) threshold 
corrections only and with the full numerical results are shown in dashed and solid lines respectively. 
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of the top-squark mixing parameters X^^ (0 and 2). All the physical masses Mj^o, M~^, 
and M3 have been set to Ms (and the same will be done for all the following plots). The dashed, 
dot-dashed and solid lines in this figure correspond to masses M/^o corrected to one-loop, two-loop 
0{as(yt) and two-loop 0{asat + a?) order[|. Fig. ^ {Xj^^ = 0) corresponds to the case of minimal 
left-right top-squark mixing, and the two-loop 0{a^) corrections are generally small, ^ 2 GeV. 
For Fig. |b (X^^ = 2), which roughly corresponds to the maximal left-right top-squark mixing 
case, we find that the two-loop Oi^a^) corrections are sizable (~ 5 GeV). 

In Fig. ^ we examine the upper limit on the Higgs boson mass M^o by including the dominant 
two-loop corrections. We show corrected massed to the two-loop O^asOtt) and 0{asCXt + a'i) orders 
in dot-dashed and solid lines for X^^ = 0,2 and the top quark pole mass Mt = 175 GeV. We see 
that maximal values for M^o of ~ 129 GeV can be reached for large tan j3 and left-right top-squark 
mixing parameter X^^ ~ 2. Without two-loop 0{af) corrections, the upper bound of M^o would 
be at ^ 124 GeV. We also show the Higgs boson masses for Mt = 180 and 170 GeV (including 
all two-loop dominant corrections) in dashed and dot-dashed lines; the masses are increased or 
decreased by ~ 5 GeV respectively. We remark that this upper bound on M^o is asymmetric with 
respect to X^^. For X^^ = —2 and Mt = 175 GeV, we find the bound is about 5 GeV lower. As 
is well know, this asymmetry arises from the two-loop 0{as(yt) corrections |P, 0. 

In Fig. 1^ we compare results from our analytical approximation formula for M^o in Sec. 2 
with those obtained by full numerical evaluations. They are shown in dot-dashed and solid lines 
respectively. The analytical approximation formula works remarkably well: it is good to a precision 
of < 0.5 GeV for almost all the parameter space. The analytical approximation has a complicated 
dependence on the yU-parameter. Numerically this dependence is quite weak: varying fi from 100 
GeV to 1 TeV for a fixed X^^ changes the Higgs boson mass by less than 1 GeV. We emphasize 
that the analytical formula is useful for several reasons: (1) the logarithmic and finite corrections 
can be easily separated, and one can weight the relative importance of these terms; (2) all terms 
can be traced back to the potential, so one can easily locate the particles giving the biggest 
contributions; (3) the formula can significantly simplify the numerical evaluations of M^o to a 
good precision. 

In Fig. 1^ we further compare the results for our RG-corrected Higgs boson masses, Eqs. ( ^8] - 
pO|), with those of the full numerical evaluation. For comparison, we have also shown two- loop 
0{asat) corrections and their RG-corrected results; they have been studied previously in 0. As 
mentioned in Sect. 3, the good agreement between these curves is an indication of the smallness 
of the logarithmic corrections beyond two-loops and illustrates the accuracy of our results. 

Finally in Figs. ^ and |^ we detail the size of the two-loop top-squark-mixing-dependent correc- 
tions in the OS-scheme and their corresponding finite threshold corrections in the RG approach. 
Fig. ^ shows in dotted lines two-loop masses without including the top-squark-mixing-dependent 
corrections of Eq. (^). Refs. @, H, have already calculated the QCD corrections, and they 
are depicted in dashed lines. The difference of the solid and dashed lines is the two-loop 0{at) 
terms which are calculated in this paper. We see clearly that these terms are sizable: for large 
mixing parameters, they increase M^o by about 4 GeV and 2 — 3 GeV for small and large tan /3 
respectively. 

Fig. 1^ shows the effect of two-loop threshold corrections A[^m\o evaluated following the RG- 
inspired analysis of Sect. 3. The dotted lines show the Higgs boson mass neglecting these correc- 

^In general, we try to follow the rule that denser lines correspond to more precise approximations. 



18 



tions; this would have been obtained by integrating two-loop RG equations (with the two-loop 
boundary threshold correction being set to zero), as we have shown in the second part of Sec. 3. 
Two-loop masses without the 0{af) threshold correction and the complete two-loop results are 
shown in dashed and solid lines respectively. The RG reshuffling of radiative corrections has ab- 
sorbed the main part of the two-loop top-squark-mixing-dependent terms of Eq. (^) into the 
RG-corrected one-loop term A[l^mf^o', the remaining genuine two-loop threshold corrections (in 
the sense of the effective field theory) are generally small, < 3 GeV. 

5 Conclusions 

In this paper we calculate radiative corrections to the lightest MSSM CP-even Higgs boson mass 
to the two-loop 0{af) order. Our analysis extends existing two-loop diagrammatic results |^, 
I, 1^ using a simpler effective potential approach and provides the most complete and accurate 
calculation presented in the literature. We also derive useful analytical approximation formulae, 
applicable when the supersymmetric particles are heavy, which accurately reproduce results from 
the full numerical study. 

Our calculation includes effects which can have an impact on the final Higgs mass but were 
neglected by previous studies. In particular, the two-loop 0{a^) top-squark-mixing-dependent 
corrections to M^o [see Eq. (pO|)] are calculated for the first term in this paper and are numerically 
important. 

We further simplify our analytical formula by reshuffling higher order logarithmic corrections 
(using RG techniques) in a compact one-loop expression [Eq. (|^)]. In that expression all mass 
parameters are evaluated at appropriate renormalization scales chosen to reproduce the numeri- 
cally most important leading and next-to-leading logarithmic corrections. The remaining two-loop 
finite terms can be interpreted as threshold corrections, and are numerically less important. This 
RG rewriting clarifies the structure of the two- loop corrections to M^o, identifies the most impor- 
tant contributions and links our work to the effective theory or RG approach, as we have shown 
in detail in Sec. 3. 

To summarize our numerical results, we have shown that two-loop top Yukawa corrections to 
M/jO are sizable for the maximal top-squark mixing case. They can increase the Higgs boson mass 
Mf^o by as much as 5 GeV (among which the top-squark-mixing-dependent corrections account 
for about 4 GeV) for small tan j3 where ht is large. The upper bound on M^o is 129 ± 5 GeV 
for Mt = 175 ± 5 GeV. Our final approximation formulae (|20| - pID and ( ^8| - |30|) have been shown 
to excellently agree with the full numerical results and can be easily implemented in precision 
numerical studies. 

Although we have focussed in this paper on the Higgs mass, it is worth mentioning that we 
have also presented in Appendix D the MSSM two-loop effective potential including top-quark 
Yukawa contributions (for general top-squark mixing parameters and any tan/3). This knowledge 
may well prove useful for other studies. 
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Appendix A: One- and two-loop scalar functions 
A.l One-loop scalar functions 

In this subsection we define the scalar functions Aq, Bq, Bi, B22 and G, which appear in one-loop 
self-energy calculations. 

The Aq function is defined by the following momentum integral in d = 4 — 2e dimensions 

,2s 



1 

m? -\- is 



m 



1 ^ . m- 



(A.1) 



where = Air/i'^e '^'^ is the renormalization scale, with 7^ the Euler constant. 
The Bq function is 

Bo(p ,m^,m2) = IGtt ji / ..^ . . -^ ^ , . .r, r:^ ^r—i 



= --/' 

e JO 



dx In 
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(A.2) 



The remaining functions can be related to and Bq as follows 
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(A.3) 
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G{p\mi,mi) = {p'-mi-ml)B,{p\miml)-A,{mi)-Ao{ml). 



(A.4) 

(A.5) 



In all one-loop expressions of radiative corrections, we adopt a (modified) minimal subtraction 
procedure to remove poles in e and keep only finite (real) parts of the above functions. 

Some useful expressions for these functions in limiting cases are (after minimal subtraction) 
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(A.6) 
(A.7) 
(A.8) 
(A.9) 
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A. 2 Two-loop scalar functions 

In this subsection we collect some useful formulae of zero-point two-loop scalar functions. They 
have been studied extensively by several groups using two different methods: a differential equation 
method |jT9|, ^ and an integral Mellin-Barnes transformation method ; their results all agree. 
Here we mainly follow Ref. [19|. 

The momentum integrals appearing in a two-loop effective potential calculation can be reduced 
to the following two types of scalar functions [corresponding to the topologies of two distinct zero- 
point two-loop irreducible Feynman diagrams (the figure-8 and sunset diagrams)]: 



J{mi,mi) = -(IGttV"'')' 



(A.IO) 



and 



(A.ll) 

The function J is symmetric in mi, m2 and / symmetric in mi,m2 and m^. 

The function J can be reduced to the product of one-loop scalar functions as 



J{ml,ml) = Ao{ml)Ao{ml) . 



(A.12) 



The function I satisfies the following first-order partial differential equation pO 
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(A. 13) 
(A. 14) 



where 

= m\ + m^ + m\ — 2m\m\ — 2m\m\ — 2m\m\ . 

This differential equation can be used to solve for the / function. The initial value of this function 
can be evaluated from ( [A.13| ) which reduces to a simple algebraic equation when m^ = mi + m2, 
i.e. R = 0. 

In our calculation, any Feynman diagram in the two-loop effective potential is subtracted by 
all its possible one-loop sub-diagrams; this is done by replacing the / and J functions as follows 
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I{ml,ml,ml) - 
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It is then straightforward to show 
J{ml,ml) = - 
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and with some effort 
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where Li2{x) is the dilogarithm function 

Lz2(a;) = 
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In the region where R^ < 0, ( [A. 18] ) should be replaced by its analytical continuation. Equivalent 
expressions for ^ also appear in and ||2^; we find that ([A.18| ) is most convenient for series 
expansions. We also define a function L for future use 

L{mj,ml,ml) = J{ml,ml) - J{ml,ml) - J{ml,ml) - {mj - mj - ml)I{ml,ml,ml) . (A.20) 

Performing a (modified) minimal subtraction (by removing the single and double poles in e), 
it is the finite (real) parts of (|A.16| ) and (|A.17D that we use in our two-loop effective potential 
expressions. We will also omit the carets of / and J to simplify the notation. 

When computing the two- loop potential, some argument of the / function, e.g. the bottom- 
quark mass mb, tree-level Higgs boson mass m/jO, can be taken to be zero. The function / is 
well-behaved in these limiting cases: 
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where we have kept only the finite terms as explained before. In ([A.21|) we have implicitly assumed 
mi > m2- The symmetry of the above expression for I{ml, ml, 0) in mi and m2 [which obviously 
follows from the definition ( A.15| ) of I] can be explicitly checked by using the identity 

^2 



Li2{x) = -Li2{x ^)--ln^( 



-X] — 



(A.23) 
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Finally, we collect expansion formulae for the function ^ which we use in the derivation of an 
analytical approximation formula for the two-loop Higgs boson mass corrections. The ^ functions 
we find can be reduced to one of the different types we list below using the relation 

^{mlmlml) = ml^{l, ml/ml ml/ml) . (A.24) 

(1) For < r < 1 and < e < 1: 

^2 



TT 



Ine - 21n(l - r) 



Inr — 2Li2(r) 



e < 2 - 21ne + Inr + 



1 + r rvr^ 



+ 



(1-r) 



^ - Ine) (1 -r^^ 



+ line - 21n(l -r) - 1 I Inr - 2Li2{r) 



r Inr 



21ne + r - 41n(l - r) 



+4rLi2 



(r) } + 0(e') 



(A.25) 



If r > 1, one uses .^(1, r, e) = r^{l,l/r,e/r) and the above expression. 
Two particular cases of the previous expansion are: 
(la) For < ei,e2 < 1: 



71 



^(l,ei,e2) = — + Inei lne2 - 2( -1 + — + Inei lne2 eie2 



71 



+ 



TT 



—2 — ^ + 2 In ei — In ei In e2 ) ei 



- Inei lei + (ei ^ e2) 



(A.26) 



with m + n = 3, and 

(lb) For < |ei|,e2 < 1: 



(81 \ 
- - -lne2je: 



+ 
+ 



2 - lne2 + (^-^ + ^\ne2 )e2 
1 1 



lne2)e?+(---lne2)et + 0(e-e^) 



(A.27) 



with m + 2n > 5. 

Finally we also give 
(2) For |ei|, |e2| < 1: 



e(l,l + ei,l + e2) = 36K +{8K- 1)6162 +[^^-^Kjele^ 

with m + n — 5. In this expansion the constant number K is given by 

K = dx ln(2cosx) ~ -0.1953256 . 



(A.28) 



(A.29) 
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Appendix B: MSSM in the leading approximation 



The general structure of the MSSM is quite comphcated, with many different fields and field 
mixings. This makes the computation of the complete potential prohibitive at two-loops. However, 
it is a good aproximation to keep only those terms of the MSSM Lagrangian which depend on 
the SU{3) gauge coupling and the top Yukawa ht (and neglect the electroweak gauge couplings 
gi , g2 and the rest of the Yukawa couplings) . We call this the leading approximation and it 
greatly simplifies our two-loop effective potential calculation. In this Appendix, we summarize 
the necessary Feynman rules for computing the two-loop potential in this leading approximation 
and also some MSSM renormalization group equations, useful to check the scale invariance of the 
potential. 



B.l Masses and Feynman rules 



The Higgs sector scalar potential in the leading approximation is 



ggs 



(m 



Hi 



(m 



Hi 



+ /i2)|iJ2|' + 5^(^1^2 + H.C. 



(B.l) 



where rriHi, rriHi ^^'^ -^m [with dimensions of (mass)^] are the soft-breaking Higgs mass parameters, 
and n the supersymmetric Higgs-boson mass parameter. Although we do not write the quartic 
Higgs couplings, which depend on the electroweak gauge coupling constants, they are responsible 
for the tree-level mass of the lightest Higgs boson, which we of course include in our calculations. 
The SU{2) doublet Higgs fields Hi and H2 can be written as follows: 



Hi 



{hi + iai)/V2 



Ho 



ht 

{h2 + ia2)/V2 



In our approximation, the mass-squared matrices for CP-even and odd Higgs fields are 



Ml 



±B„ 



±B„ 



m 



H2 



(B.2) 



(B.3) 



where the positive (negative) sign applies to the CP-even (odd) fields respectively. The charged 
Higgs fields have the same mass-squared matrix A^^ as the CP-odd Higgses. 

The CP-even interaction eigenstates hi, h2 are rotated by the angle a into the mass eigenstates 
H^ and h^. Similarly, the CP-odd states ai, 02 (charged states hf, h^) are rotated into mass 
eigenstates and {G~^ and H^) by the angle (3. This angle f3 is conventionally defined in terms 
of the CP-even Higgs field VEVs, (/ii,2) = "^1,2, by tan/3 = V2/V1. The fact that P diagonalizes 
is obvious when the minimization conditions of the potential ( [B.l|) , +11'^ = — i?^tan/3 
and m|^2 + = —B^ cot /?, are imposed and the soft parameters in the matrix are replaced by 
tan j3 and m\o = —B^(ta.n [3 + cot /?). Since we have neglected all gi, g2 related terms in ( p.l|) , in 
our approximation (we use shorthand notations cp = cos/5, = sin/3, etc.) 



3/3, 



and Sr 



(B.4) 



This approximation is excellent when M^o 3> Mz but would fail for M^o ~ Mz- The effect is 
numerically relevant for the tree level masses and we take it into account, but it may be consistently 
neglected in the two-loop corrections. 
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The (field-dependent) top and bottom squark mass-squared matrices (neglecting the D-terms) 

eg 



Mi 



Mt 



-^ht{Ath2 + fihi] 
Ml + \hlhl 



-^ht{Ath2 + ^ihi 
Ml + \hlhl 

■^hb{Ahhi + /i/i2 



Ml + \hlhl 



(B.5) 



(B.6) 



-^hb{Abhi +/i/i2) 

where M^, M^ (M^) are soft-breaking mass parameters of the left- and right-handed top(bottom)- 

squarks Q and U (D); At and Ab are the usual trilinear soft-breaking parameters. We denote the 
mass eigenvalues of the matrix ( |B.5| ) by m^^ , m^^ the mixing angle by 6^, and the corresponding 
quantities for the matrix ( |B.6| ) by m^^, and 6'^. 

The Feynman rules for Higgs/Goldstone-boson-squark trilinear coupling are simply —iX, with 
A as listed below: 



A, 



A 



H+i2b 



A 



H+t2b2 



htCfB 

htCp 



{ctmt + StYt)cb + rribSbSt 
{ctmt + StYt)sb - rUbCbSt 
{stnit - CtYt)cb - nibSbCt 
{stnit - CtYt)sb + rribCbCt 



-hbSp 
+hbSp 
+hbSp 
-hbSp 



{cbrrib + SbYb)ct + rritStSb 
{sbTTib - CbYb)ct - mtStCb 
{cbnib + SbYb)st - mtCtSb 
{sbTTib - CbYb)st + mtCfCb 



(B.7) 



and 



-^G+tifci = -htSfsiCtTTlt + StXt)Cb + hbCfsiCbTTlb + SbXb)Ct , 

^G+hh = htSp{ctmt + StXt)sb + hbCp{-Sbmb + CbXb)ct , 

-^G+tafci = htSp{stmt - CtXt)cb - hbCp{cbmb + SbXb)st , 

^G+i2b2 = htSp{-Stmt + CtXt)sb + hbCp{sbmb - CbXb)st , (B.8) 



^HH^h = ^y^htimt + StCtY")sa , >^H0t2t2 = \/2/it(mt - 
Kohh = V2ht{mt + StCtX")ca , Kohh = "/^htirrit - StCtX")c, 



a 5 
a ; 



^H0hi2 — —7^htC2tY"Sa , \oiJ^ — —7=htC2tX^Ca 



■^AOtita — ~^m2h ~ ~^^tYtCf3 ) -^GOtita — --^GOtaii " —^hfXtS/s , (B.9) 



where q = cos 6^, St = sin 6^, C2t = cos 26^ (with similar expressions for 6^ functions) and 

Xt = At + fi cot p, Yt = At-fitanp, (B.IO) 

Xb = Afe + /itan/?, Yb = Af,-/icot/3. (B.ll) 

^In this revised version, we have also included bottom Yukawa terms in the Feynman rules, they will be used in 
the expanded two-loop effective potential expression ( D.6 ). 
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In addition, we find convenient to define the a-dependent quantities 



X^" = At-fitana, F^" = At + fi cot a, (B.12) 

= Ab — ficota, = Ab + fitana, (B.13) 

which tend to the corresponding quantities without the a label (X" Xt, etc) in the limit 
TTlA > Mz. 

Couplings similar to the above ones but for bottom squarks can be obtained directly from 
(|B.7|) , (|B.8|) and (|B.9|) : simply make everywhere the replacements {ht ^ hb, rrit ^ rUb, 9i ^ 
01, Xi"^ ^ Xi"\ f/"^ ^ f/"^} and {c„ ^ ^ -s^} for the couplings to {H+,A^,H^} or 

{ca ^ —Sa, Cj3 ^ S/3} for the couplings to {G~^, (7°, /i"}. 

The couplings of squarks to neutralinos and charginos are very simple in the leading approxi- 
mation, since the gaugino-Higgsino mixing can be neglected and the only interactions are Higgsino- 
squark interactions. The Feynman rules for the h^ttj couplings can be written as —i{aVL + bV^) 
and that of h'^thi as iC~^{aVL + bVu) {Vl,r are chiral projectors and C the charge-conjugation 
matrix), with 

_ • -h - -h _ ^* 
- 7, - - 

~'^h\tt2 - ^^hltt2 - ^h\tix - '^"hOtii - -^^t ' 

H+tbi^ = , a-^+f,^^ = -htst , a^+^^-^ = -htCt , (B.14) 

when /i > 0; for /i < 0, we only need to interchange a^o^^. and a/jOf^., as well as 6^0^^. and b^^^,. 

Other Feynman rules of Olg^) and 0{ht) vertices are exactly the same as in the general MSSM 
and we do not present them explicitly. 



B.2 Renormalization group equations 



The MSSM RGBs |23| that we will use to check the invariance of the potential to two-loop order 
under renormalization scale transformations are the following. First, we need the two- loop RGEs 
for those parameters entering in the tree-level potential ( p.l|) 

^ - -^M\ + ;\' {M\ + 2M| - 2M^At) - . ^ ' {M\ + A\) , (B.15) 



4 



(91n/i _ d\nh2 _ 3h^t Sg^h^t 9 



dlnQ^ dlnQ^ 327t^ {IGTT^y 2 (167r2)2 



(B.16) 



ainQ2 I67r2 V 2 '^7 (167r2)2 V 2 

^^'^+2^/.) , (B.17) 



(167r2)2 V 2 

where Aii= fnjj^ + M~ + M~ + Af. Then we need one-loop RGEs for those masses entering in 
the one-loop potential 
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167r' 



+ K 



'im\ + (1 + s^^)M\ + Zs^mt ^At + ^ j 



167r' 



dm 



t2 



Xt 



+ hi 



3ml + (1 + cl)Ml - 3s2tmt At + 



3Xt 



IGtt 



1677 



2 ^"^HO 

dlnQ^ 



3h 



+ Atii 



(B.19) 

(B.20) 
(B.21) 
(B.22) 



where = s^, c^, s^, and -E„ = S2a, —S2a, — 52/3, S2/3 for n = 1, 2, 3, 4. [Here we use the a angle 
to keep track the and contributions; it can be replaced by the (3 angle as in ( p.4|) in the 
leading approximation.] The ordering of the Higgs/Goldstone bosons are = H^,h^,G^ and 
A^ for n = 1,2,3,4 and H+ = G+,H+ for n = 1,2. Eqs. (p:T9|^:2^ ) seem unfamiliar, but they 
follow directly from ( [B.3| ), ( |B.5| ) and the one-loop MSSM RGEs of soft parameters entering those 
equations. 

Using ( |B.19| ) and ( |B.20D , we find one-loop RGEs for Xt and m?, the (arithmetic) average of 
the (squared) top squark masses. They are 



IGtt 



dXt 
(91nQ2 

(91nQ2 



16 

-glMs + 3hl{At + Xt) . 
-fgliml + Mi) + hl(3ml + lMl, , 



(B.23) 
(B.24) 



these two equations are used in Sec. 3 for the RG discussion of the formula for the Higgs boson 
mass M/jO. Eq. ( |B.23| ) can also be derived from ( |B.1(]| ) and one-loop RGEs of At, fi and tan/5. 



Appendix C: One-loop self-energies 



In this appendix, we collect formulae for those MSSM one-loop self-energies which are necessary 
for our analysis. We present these self-energies in the leading approximation of keeping only ht 



and ^fs-dependent terms, as explained in Appendix B; their full form can be found in Ref. jlO 
which we follow for notation. (See also for top quark/squark self-energies.) 
-Top quark: 



167r2 Et(p2 



4^ 
3 

S2t m-g 



hi i 2 



mt B,{p\mlml) + B,{p\mlml) 
5o(/, m\) - Bo{p^, m?, m| 
2Bi {p^ , ml , m% ) + Bi {p^ , ml , m\o 



m; 



-mi(5-31n^ 
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+ St 



2B,{p\mf,mi.) + Bi{p\mlm'z) 



+ B, /i^, ml ) + 5i m| ) + Bi /i^, m?^ ] 



(C.l) 



where we have assumed all heavy Higgs bosons have mass m^o much larger than the masses of 
the light Higgs and W^-boson, taken to be ~ mz- 

From ( (J.l ) we find the running top-quark mass at the scale Q (under the simplified assumptions 
of a common heavy SUSY scale Ms while the n parameter is left free, see Sec. 2) 



mUQ) 



(C.2) 



+ 



3hl 
327r2 



(1 + 4) --In -f +4 --In -i 



Q 



l-/i2 



1 — fi'' 



In this equation we have neglected the external momentum and used ( |A.3|) and ( |A.6| ). We have 
used the reduced parameters Xt = Xt/Ms, fi = fJ^/Ms and Mf is the top quark pole mass (we use 
capital letters to denote on-shell mass parameters). 

-Top squarks: 



M|G(p2,^,^2)^^2 



+ s 



A,{ml)-{p' + ml)Bo{p\ml,0) 
Ao(m|)-(/ + m|)5o(/,m|,0) 



+ h' 



n=l 



4 2 



n=l 1=1 



n=l 



16n' Urrip') 



Ao(m?)-(p2 + m?)i?o(p^m? 0) 



+ ct 



+ h' 



Aoiml) - ip' + ml)Bo{p',ml,0) 
CtMK) + StAiml) + Ao(m?J 



+ -J2DnAo{mlo) + J2Dn+2Ao{ml+) + G{p'',ij\m'^)+G{p^,fi\ml) 

^ n=l " n=l 



4 2 



71 = 1 1 = 1 



n=l 



S2t{v' + ml)Boip\ ml , 0) + S2t(p' + m?JEo(p^ m|, 0) 



t2> 
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+ 4mgmtBo{p ,mg,m^ 

4 2 

n=l i=l 
2 



Ao(m|) - Ao{m 



(C.5) 



n=l 



where XuH^h = Ct\HOhh - ^t^i?o?2?i' ^HOt^ti = sAH^i^h + (^t^ROhh^ etc.. The symbols are 
defined after ( |B:T0| ). 

From (|C.3[ - |C.5| ) we derive relations between running and on-shell top-squark masses and mixing 
parameters using the following one-loop relationships (for = = 1/2): 



- M^ + mi + m,X, - -Re [n,-^,-^ (M/J + H,-^,-^ (M^ ) J - Re H,-^,-^ (M^ ) , 



+ - m^Xt - -Re 



%f,(M|) + n,^,jM|)J + Re %,,(M|) , (C.6) 
we obtain (assuming again a common heavy SUSY scale Ms and leaving free the //-parameter): 



m}{Q) 



Mil 1 - 



+ 4 1 



-r^-in^^ 



37r2 

TT 



+ 



327r2 



Mi 



{Xts$ + Y,^c$) (2-ln^ 



+ /i^ln/i^ + (1 -/i^) 3_21n— f - {I - j^y - jl 



mtXt{Q) 



M'^ ^ A/f2 

4(2-ln|f)+2X,ln|f 



(C.7) 



--.2,^/. 2-ln 



i^t4 + xJi 



M| 



l-jl^ + fi^\nfi^ + {l-fi'^) ln(l - /i^ 



+ 



7r \ 



2 3v^/ ^ * 2 ^ * * V ^5 / 3 



2 i>2 ■ 



s^X/X( In 2 



where we have used ( |A.7| ), ( |A.9| ) and the definition Yt = (A — fita.n P)/Ms. 
-W boson: 



B,,{p^mlml)-lMml 



2si 



B,,ip^mlml)-tMml 



(C.9) 



This gives (under the assumption of the simplified SUSY spectrum of Sec. 2, described already 
for previous self-energies) 



4(g) = -[M2, + Ren^^(M2,)] 



2„2 



1 - 



327r2 



mr 



-61n^ + 3 + X/ 



(C.IO) 
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where we have neglected the external momentum in ( |C.9| ) and used ( [A.3[ - [A.6D . 

-Higgs boson: We need only the difference 



167r^ 



o 7,2 2 2 



Bo{0, ml, ml) - Ami ^^ob^ ml, ml] 



p2=0-' 



d 



(C.ll) 



p2=0 



where A/jj.t. are defined in ( [B.9|) . Using ( |A.6D and (|A.7| ) we get (|T9l). 



Appendix D: MSSM effective potential to the two- loop or- 
der 

In this Appendix we present the MSSM effective potential for the (real) neutral components of the 
Higgs fields up to the two-loop level in the leading approximation (which neglects all dimensionless 
couplings except ht and g^). We first write the potential as 

V{h, h^) = Kac + ^(/il, + V^{h^, h^) + F2(/li, h^) , (D.l) 

where V^ac is a field-independent vacuum energy termj]. The tree-level potential Vq is 

H(/ii, h2) = \{ml^ + ^?)hl + i(m|^ + ^?)hl + B^h^h2 , (D.2) 

which simply follows from substituting Eq. ( p.2| ) into the MSSM Higgs sector scalar potential 

The one-loop potential is well known and the 0{asOtt) part of the two-loop potential was 
computed in we list them here for completeness and for future reference. The complete one- 
loop potential in Laudau gauge isQ 

lQT^^V,{h,M) = Y^^cIT. H{ml)~2Him})\+3H{m'^) + ^H{ml) 

+ lj:m<o) + j:H{ml,)-2j:H{ml,)-j:H{mlo), (D.3) 

n=l n=l i=l 1=1 

where / sums over all the (s) quarks and (s)leptons, is the color factor, 3 for (s) quarks and 
1 for (s)leptons. Following the leading approximation, we keep only the numerically important 
parts, i.e., those from top (s)quarks. In Eq. ( |D.3| ), xt{i = 1,2) and ^^(^ = 1)2,3,4) represent 
charginos and neutralinos, and the function H is 

rr/ 2\ m^ (^ 3\ ,^ ^. 

^("^ ) = ^ In— - - . (D.4) 



Q 



^This term is a function of the soft-breaking parameters; it is needed for the invariance of the potential under 
a RG transformation. 

8 We adopt the (modified) DR -scheme of Ref. ||. 
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The QCD contribution to the two-loop effective potential in the MSSM is 



+ st) E JK, m|) + 2slcp{ml,ml) + m|/(m| , m| , 0) 



i=l 



+ ^L{m\,rnl,ml) - Am-g mt StCt I{m\,ml,ml) - I{m\,rnl,ml 



i=l 



+ 



(D.5) 



where g is the gluino, with tree-level mass given by the SU{3) gaugino soft mass M3. The last 
term, obtained by interchanging variables, gives the contribution from sbottoms. Note that there 
is no mixed contribution involving stops and sbottoms, even if such mixed couplings exists [from 
SU{3) D-terms]. The two-loop scalar functions J, J and L in Eq. ( p.5| ) are given in Appendix A, 
Eqs. ( [ATtD , ([AT6D and (^:20|).P| 
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Figure 8: Feynman diagrams for the two-loop effective potential of order 0{a'^) in the MSSM. 
represent if", h^, and for n = 1, 2, 3, 4. represent G"^, if"*" for n = 1,2. The neutral and 
charged Higgsinos (= X34) ^-iid (= x^) have degenerate mass of 

The top Yukawa contribution to the two-loop potential is a new result of this paper. The 
relevant Feynman diagrams are shown in Fig. ^. To simplify the final result, we neglect left-right 
mixings in the bottom-squark sector and the gaugino-Higgsino mixings in the neutralino-chargino 
sector (under this assumption, the Higgsino masses are simply these simplifications are valid 
in the leading approximation. Using the Feynman rules given in Appendix B, we find (the last 

^The procedure we have followed of subtracting all possible one-loop subdivergences to define these functions 
is an alternative to the direct way used in Ref. The direct way requires the computation of some one- loop 
quantities to order 0{e); perhaps we find the subtraction method simpler. We have explicitly checked that, in the 
particular limit studied in we exactly reproduce their unexpanded mass formula, Eq. (11) of which shows 
the equivalence of both methods. 
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diagram of Fig. ^ is of order but does not contribute to m/jo) the top and bottom Yukawa 
contribution to the two-loop potential^°: 

{IGny {V2t{hi,h2) + V2b{h,h2)) = 

-IT L{rn]j^^,ml, ml) ± 2m^J{mj^o,m'f, ml) + Y,J (m|, m^o; 



.n=l 



i=l 



+ sl J (m| , m\+ ) + cp (m| , + ) + J (m?^ , ) + J (m?^ , ) 



n=l 



+3|(/i?s2 + hlcl)L{ml,fi^ml) + {h^cl + hlsl)L{ml, ^\ml) 

-2fimbhthbS2t I{mf^,fi'^,ml)-I{mf^,n'^,ml) | 

+ {/it ^ rrii ^ mb, m^^ ^ m^^, ^ ^ D^}] 

+3|(/it 4 + hlcl)L{ml+,ml, ml) + {h]cl + hlsl)L{mjj+,ml, ml) 



+2mtmbhthbS2f3 



2 4 



I{ma+,mt,mt^) - I{mjj+,mt,mf,) 

2 2 



(D.6) 



i,j=l n=l 



i,j=l n=l 



where, in the first hne of Eq. ( |D.6| ), positive and negative signs apply to CP-even and odd 

{A^, G'^) Higgs/Goldstone bosons respectively, and in the last hne fi = {U,bi}. We also have, for 
n = 1,2, 3, 4, Dn = {s^, c^, s|, c|} and D'^ = {c^, s^, c|, s|}. The ordering of the Higgs/Goldstone 
bosons is = H°, h°, G° and A° for n = 1, 2, 3, 4 and i/+ = G+, H+ for n = 1, 2. 

Two tests can be applied to check the correctness of the effective potential V{hi,h2). First, 
the potential should vanish in the supersymmetric limit [i.e., when all soft-breaking parameters 
are taken to be zero), and second, the potential V{hi,h2) should be invariant under changes of 
the renormalization scale, up to the order of our perturbative calculation. The vanishing of the 
potential in the supersymmetric limit is proved by simple algebra. In the following we show the 
invariance of the two-loop potential under a RG transformation. 

Using the derivatives of /, J and L functions with respect to the renormalization scale Q 



dl{mi, m^, ml 



91nQ2 

i) 



Ao{mf)+ml 



dJ{ml, m^ 



(91nQ2 



= -E 

i=l 

m\AQ{m\) + m\AQ{m\) 



(D.7) 



:d.8) 



In this revised version, we have also included the bottom Yukawa contributions for completeness. All analyses 
in the main text use only the top Yukawa contributions as in the previous version. 
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dL{m\^ 1712, ^\ 



+ 771, 



2\2 



(7712 + 7713) 



:d.9) 



and the one-loop MSSM RGEs for top-(s)quark and Higgs boson masses, Eqs. ( [B.18| - p^ ), we 
find that the RG variation 



dV2 
9InQ^ 



(IGtt , 
fl67r2) 



(D.IO) 



modulo terms independent of the Higgs field h^- Here V^^^Vi represents the one-loop RG variation 
of the one-loop potential Eq. ( p.3|) . This result agrees exactly with the two-loop RG variation of 
the tree-level potential [cf. Eqs. (^) and ( P7[3| - |BT^ )], so that 



dlnQ 



91nQ2 



. 



(D.ll) 



Note that this is a nontrivial check that all In terms cancel with each other between Eq. ( p.lO| ) 
and "D^^^Vo; this cancellation guarantees the correct leading and next-to-leading order logarithmic 
behavior of the effective potential. 

To derive the analytical expression of Atti^o in Sec. 2, we need to expand the two- loop potential 
V2 in powers of rrit/Ms and rritXt/M^; besides many straightforward expansions, we have used 
(|A.25|) - (|A.28| ) for the t — q — h and t — q — h diagrams of (|D.6| ), with q = t or b. 
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